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GATED CONE-BEAM COMPUTED TOMOGRAPHY RECONSTRUCTION WITH MOTION COMPENSATION 



5 The present invention relates to the field of computed tomography such 

as cone-beam CT. In particular, the present invention relates to a method of 
reconstructing projection data from a gated projection data set, to an image processing 
device, to a computed tomography apparatus and to a computer program for 
reconstructing projection data from a gated projection data set. 

10 

Cone-beam computed tomography (CBCT) enables volumetric imaging 
at high spatial resolution. With the improved acquisition speed, cardiac CT imaging has 
become feasiible. However, due to the gating of the projection data, the source orbit is 

15 interrupted, resulting in a reduced dose efficiency and prohibiting the use of exact 
reconstruction techniques as, for example, described in B.D. Smith "Image 
reconstruction from cone-beam projections: necessary and sufficient conditions and 
reconstruction methods" IEEE trans, med. image., MI-4: 14-25, 1985. 

In such cardiac cone-beam computed tomography systems, the source 

20 path i.e. the path of the source during which data is gathered is interrupted due to the 
fact that the projection data is gated. In cardiac cone-beam computed tomography, the 
gating may be performed in accordance with an electro-cardiogram (ECG) or any other 
suitable means for detecting the movement of the heart. Due to the gating, a majority of 
projection data is discarded which strongly decreases the dose utilization. Furthermore, 

25 in some cases, the 3D completeness criterion for exact reconstmction as formulated by 
Tuy may be violated resulting in incomplete radon data and thus prohibiting the use of 
exact reconstmction techniques. 

30 It is an object of the present invention to provide for an improved 

reconstruction. 
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According to an exemplary embodiment of the present invention as set 
forth in claim 1, the above object may be solved by a method of reconstructing 
projection data from a gated projection data set wherein the gated projection data set is 
firstly acquired. The source trajectory, i.e. the trajectory of the source of radiation used 
5 for acquiring the gated projection data set, has at least one gap due to gating. In case the 
method is applied to a CBCT, the gating is preferably performed in accordance with a 
movement of the heart which, as described above, may be detemiined by means of an 
electro-cardiogram or other suitable means. Then, new projection data is determined 
corresponding to the at least one gap and the gated projection data is supplemented with 

10 the new projection data to compensate for the at least one gap in the source trajectory. 

Advantageously, the supplemented gated projection data may thus be 
completed in a way that, for example, the three-dimensional completeness criterion for 
exact reconstruction is no longer violated resulting in complete radon data. This 
advantageously may then allow the use of exact reconstruction techniques. 

15 According to another exemplary embodiment of the present invention as 

set forth in claim 2, a four-dimensional image data set is reconstructed from the gated 
projection data by using a cone-beam computed tomography reconstmction. On the 
basis of this four-dimensional image data set as set forth in the exemplary embodiment 
of claim 3, a four-dimensional vector field may be determined which describes the 

20 motion of the object of interest which is imaged. This motion field may then be used for 
performing a motion compensation of the gated projection data. 

In other words, an approximate reconstruction technique may be used to 
reconstruct a four-dimensional data set. From this four-dimensional data set, a four- 
dimensional vector field is calculated. Then, the four-dimensional image data is motion 

25 compensated. 

According to another exemplary embodiment of the present invention as 
set forth in claim 4, the four-dimensional image data is motion compensated and then, 
the motion compensated four-dimensional data is used to calculate new projection data 
for completing the original data set, i.e. for filling the gap. Advantageously, this may 
30 allow to form a complete gated projection data set without gaps which may allow for 
the application of exact or quasi-exact reconstmction techniques. 

Then, as set forth in the exemplary embodiment of the present invention 
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as defined in claim 5, an approximate or an exact reconstruction algorithm or method 
may then be applied to the gated projection data which was supplemented with the new 
projection data which allows to reconstruct the final image data without interruption for 
the image generation. 

5 Advantageously, this may allow for an improved image quality and for a 

fast and robust reconstruction. 

According to another exemplary embodiment of the present invention as 
set forth in claim 6, an image processing device is provided comprising a processor for 
reconstruction projection data from a gated projection data set. Advantageously, the 
10 image processing device is adapted to perform an operation to supplement the gated 
projection data such that gaps in the data caused by an incomplete source trajectory are 
filled. 

Advantageously, this image processing device allows for a very exact 
and fast reconstruction. Another exemplary embodiment of the image processing device 
15 is set forth in claim 7. 

According to another exemplary embodiment of the present invention as 
set forth in claim 8, a computed tomography apparatus is provided comprising a 
processor which is adapted to perform an operation in accordance with the method of 
the present invention. Advantageously, this computer tomography apparatus allows, for 
20 an exact or quasi-exact reconstruction of data in spite of an interrupted source trajectory 
due to a gating. 

According to another exemplary embodiment of the present invention as 
set forth in claim 9, a computer program for a data processor for reconstructing 
projection data firom a gated projection data set is provided. The computer program 

25 according to the present invention is preferably loaded into a working memory of a data 
processor. The processor is thus equipped to carry out an exemplary embodiment of a 
method of the present invention. The computer program may be stored on a computer 
readable medium, such as a CD-ROM. The computer program may also be presented 
over a network such as the Worldwide Web, and can be downloaded into the working 

30 memory of a data processor firom such a network. The computer program may be 
written in any suitable programming language, such as C-H-. 

It may be seen as a gist for an exemplary embodiment of the present 
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invention that projection data may be obtained corresponding to an uninterrupted 
source trajectory from a gated data set where the source path was interrupted due to the 
fact that the data acquisition was gated. For this purpose, a motion compensation is 
applied. An approximate reconstruction is used to reconstruct a four-dimensional data 
5 set. From this data set, a four-dimensional vector field is calculated. The four- 
dimensional image data is then motion compensated and subsequently used to calculate 
new projection data completing the data set. The new projection data corresponds to the 
missing data, i.e. to the gaps in the interrupted source path. Then, an appropriate or 
even exact reconstruction algorithra may be applied to the data set without interruption 
1 0 for the image generation. 

These and other aspects of the present invention will become apparent 
from and will be elucidated with reference to the embodiments described hereinafter. 

Exemplary embodiments of the present invention will be described in the 
following with reference to the following drawings: 

15 

Fig. 1 shows a schematic representation of an exemplary embodiment of 
a cone-beam computed tomography scanner as it may be used for cardiac 
cone-beam CT according to the present invention. 
20 Fig. 2 shows an exemplary embodiment of a method of operating the 

computed tomography apparatus of Fig. 1. 

Fig. 3 shows a schematic representation for further explaining a principle 
of an exemplary embodiment of the present invention. 
Fig. 4 shows a simplified schematic representation of an exemplary 
25 embodiment of a data processing device according to the present 

invention. 

Fig. 1 shows an exemplary embodiment of a computed tomography 
30 apparatus according to the present invention. Preferably, this is a cone-beam computed 
tomography apparatus (CBCT) where a cone-beam 6 is applied to an object of interest. 
With reference to this exemplary embodiment, the present invention will be described 
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for the application in cardiac cone-beam CT. However, it should be noted that the 
present invention is not limited to a cardiac CT, but may be applied to any CT imaging 
method where an incomplete data set is used due to a gating or a suitable sampling. In 
other words, the present invention may be applied to the imaging of moving objects in 
5 general where certain movement stages are imaged. 

The CBCT scanner depicted in Fig. 1 comprises a gantry 1 which is 
rotatable around a rotational axis 2. The gantry 1 is driven by means of a motor 3. 
Reference character 4 designates a source of radiation, such as an x-ray source which 
according to an aspect of the present invention emits a polychromatic radiation. 
10 Reference character 5 designates a first aperture system which forms the 

radiation emitted from the radiation source 4 to a cone-shaped radiation beam 6. 

The cone-shaped radiation beam 6 (or cone-beam) is directed such that it 
penetrates the object of interest 7 arranged in the centre of the gantry, i.e. in an 
examination region of the CBCT scarmer, and impinges onto the detector 8. As may be 
15 taken from Fig. 1, the detector 8 is arranged on the gantry 1 opposite to the source of 
radiation 4. The detector 8 depicted in Fig. 1 has a plurality of detector lines each 
comprising a plurality of detector elements. 

The detector lines of the detector 8 are arranged at the gantry 1 such that 
the lines are perpendicular to the rotational axis 2. Furthermore, columns of the detector 
20 8 are essentially parallel to the rotational axis 2. In other words, the detector 8 may be a 
two-dimensional detector. 

The apertures of the aperture system 5 is adapted to the dimensions of 
the detector 8 such that the scanned area of the item of interest 7 is within the cone- 
beam 6 and that the detector 8 covers the complete scanning area. Advantageously, this 
25 allows to avoid unnecessary excess radiation applied to the object of interest 7. During 
a scan of the object of interest 7, the source of radiation 4, the aperture system 5 and the 
detector 8 are rotated along the gantry 1 in the direction indicated with arrow 16. For 
rotation of the gantry 1 with the source of radiation, the aperture system 5 and the 
detector 15, the motor 3 is connected to a motor control unit 17 which is connected to a 
3 0 calculation unit 1 8 . 

In Fig. 1, the object of interest 7 is disposed on a rest 19, which may be 
movable. According to an exemplary embodiment of the present invention, as 
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mentioned above, a circular data acquisition is performed where the x-ray source 4 is 
displaced along a circular source trajectory i.e. is rotated in a rotational plane around 
the rotational axis 2 without a movement of the object of interest in a direction parallel 
to the rotational axis 2. Thus, the object of interest may be immobile when the x-ray 
5 source 4 performs a circular movement. However, according to another exemplary 
embodiment of the present invention, for achieving a helical source trajectory, the 
object of interest 7 may be moved on or together with the rest 19 in a direction parallel 
to the rotational axis 2. 

The detector 8 is connected to a calcialation unit 18. The calculation imit 

10 18 receives the detection results i.e. read-outs from the detector elements of the detector 
8 and determines a scanning result on the basis of ttie scanning results from the detector 
8. In addition to that, the calculation unit 18 commixnicates with the motor control unit 
17 in order to coordinate the movement of the gantry 1 with the motor 3 and 20 or with 
the rest 19. The calculation unit 18 is adapted for reconstructing an image from read- 

15 outs of the detector 8. The image generated by the calculation unit 18 may be output to 
a display (not shown in Fig. 1) via an interface 22. 

Fig. 2 shows an exemplary embodiment of a method of operating the 
computed tomography apparatus of Fig. 1 . 

After the start in step SI, the method continues to step S2 where a data 

20 acquisition is performed. In particular, in step S2 projection data of a three-dimensional 
volume is acquired. Due to a gating process perfonned, the three-dimensional volume 
that is acquired has gaps. For example, for a cardiac cone-beam CT, the source path is 
intermpted due to the fact that the projection data is gated according to an 
electrocardiogram (ECG) or to any other suitable means adapted to determine a 

25 movement of the heart. Due to this, a part of the projection data that is actually acquired 
is discarded which strongly decreases the dose utilization. This causes the gaps in the 
source trajectory i.e. the gap in the three-dimensional image volume that is acquired. 
Furthermore, due to this, a 3D completeness criterion for extract reconstruction may be 
violated which results in an incomplete radon data set. 

30 Then, in the subsequent step S3, a four-dimensional image data set is 

reconstructed. For the case of a cardiac CBCT reconstruction, a four-dimensional image 
data set may be reconstructed by means of an approximate reconstmction algorithm 
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such as, for example, described in Grass, M. et. aL "Helical cardiac cone-beam 
reconstruction using retrospective ECG gating" Phys. Med. Biol. 2003, which is hereby 
incorporated by reference or KachelrieB, M. et. al. "ECG-correlated image 
reconstruction from subsecond multi-slice spiral CT scans of the heart" Med. Phys., 
5 27(8): 1881-1902, 2000, which is also hereby incorporated by reference. 

Then, in a subsequent step S4, a vector field is calculated. In particular, 
from the four-dimensional image data reconstructed in step S3, a four-dimensional 
vector field is calculated using, for example, a three-dimensional registration technique 
or block matching algorithm, such as described in Schaffter T. et. al. "Motion 

10 compensated projection reconstruction" 41: 954-963, 1999, which is hereby 
incorporated by reference. 

Then, in subsequent step S5, the four-dimensional vector field 
determined in step S4 is used for a motion compensation of the three-dimensional 
volume acquired in step S2. Then, in the subsequent step S6, the motion compensated 

15 image data is used to fill gaps in the trajectory by calculating new projection data with 
e.g. forward projection. 

Then, in subsequent step S7, an approximate or an exact reconstruction 
algorithm, such as described, for example, in Katsevich, A. "Analysis of an exact 
inversion algorithm for spiral cone-beam CT" Phys. Med. Biol., 47: 2583-2597, 2002, 

20 and Katsevich, A. "Theoretically exact FBP-type inversion algorithm for spiral CT" 
SIAM J. App. Math., 62: 2012-2026, 2002, which are both hereby incorporated by 
reference, may be applied. Such reconstruction algorithm may be used to reconstruct 
the final volume (images) using the motion compensated projection data. 
Advantageously, this may allow for an improved dose efficiency and thereby a decrease 

25 of a dose of radiation applied to, for example, a patient. 

Advantageously, the above method may ftirthermore allow for the 
application of large area detectors for which an exact or quasi-exact reconstruction 
approach, such as the one described above, is necessary. 

Fig, 3 shows a simplified schematic representation for visualizing the 

30 method described with reference to Fig. 2. As may be taken from Fig. 3, due to the 
gating, the source trajectory from which the data is acquired which is used for 
subsequent reconstruction is interrupted. In the case of cardiac CT, the source trajectory 
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is interrupted due to ECG-gating. According to the present invention, in the subsequent 
step 83, an approximate three-dimensional reconstruction is performed for different 
phases. Then, the method continues to steps S4 and S5 where a four-dimensional vector 
field determined from four-dimensional image data is used for motion compensation of 
5 the image data. Then, the motion compensated data is used to fill the gaps in the 

trajectory by calculating new projection data with forward projection. Then, the motion 
compensated and supplemented image data is used for an approximate or an exact 
reconstruction algorithm which allows, for example, for an exact three-dimensional 
reconstruction from the completed data set. 

10 Fig. 4 shows an exemplary embodiment of a data processing device, such 

as an image processing device, for performing the method described with reference to 
Figs. 2 and 3. As may be taken fi:om Fig. 4, a central processing xmit (CPU) or image 
processor 51 is connected to a memory 52 for storing a gated projection data set, any 
intermediate data or the finally reconstructed data. The data may be acquired by a 

15 CBCT scanner such as the one depicted in Fig. 1 . For this, the image processor 5 1 may 
be connected to such a CBCT scanner and/or to a plurality of input/output/network or 
other diagnosis devices. The processor 1 is furthermore comiected to a display 54 (for 
example to a computer monitor) for displaying information or images computed or 
adapted in the image processor 1 . An operator may interact with the data processor 5 1 

20 via a keyboard 55 and/or other input or output devices which are not depicted in Fig. 1 . 

The present invention described above may, for example, be applied in 
the field of medical imaging. However, as described above, the present invention may 
also applied in other areas where moving objects are to be examined, such as in the 
field of non-destructive testing. 



